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Abstract. In this work we derive a lower bounds for the Hausdorff and fractal dimen- 
sions of the global attractor of the Sabra shell model of turbulence in different regimes of 
parameters. We show that for a particular choice of the forcing and for sufficiently small 
viscosity term v. the Sabra shell model has a global attractor of large Hausdorff and fractal 
dimensions proportional to log^^v"' for all values of the governing parameter e, except for 
E = 1 . The obtained lower bounds are sharp, matching the upper bounds for the dimension 
of the global attractor obtained in our previous work. Moreover, we show different scenar- 
ios of the transition to chaos for different parameters regime and for specific forcing. In 
the "three-dimensional" regime of parameters this scenario changes when the parameter e 
becomes sufficiently close to or to 1 . We also show that in the "two-dimensional" regime 
of parameters for a certain non-zero forcing term the long-time dynamics of the model 
becomes trivial for any value of the viscosity. 
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1. Introduction 

Shell models of turbulence have attracted interest as useful phenomenological models 
that retain certain features of the Navier-Stokes equations (NSE). In this work we continue 
our analytical study of the Sabra shell model of turbulence that was introduced in 1141 . For 
other shell models see 0, ||9l, HI, 1171 . A recentreview of the subject emphasizing the ap- 
plications of the shell models to the study of the energy-cascade mechanism in turbulence 
can be found in J2|. 

The Sabra shell model of turbulence describes the evolution of complex Fourier-like 
components of a scalar velocity field denoted by u„. The associated one-dimensional 
wavenumbers are denoted by k„, where the discrete index « is referred to as the "shell 
index". The equations of motion of the Sabra shell model of turbulence have the following 
form 

^ i{ak„+iu„+2U*+i+bk„u„+iu*_i -c^„_im„_im„_2) -V^,^m„ +/„, (1) 

for « = 1,2,3, .. ., with the boundary conditions m_i = uq = 0. The wave numbers k„ are 
taken to be 

k„=koX'\ (2) 

with A- > 1 being the shell spacing parameter, and ko > 0. Although the equation does not 
capture any geometry, we will consider L = k^^ as a fixed typical length scale of the model. 
In an analogy with the Navier-Stokes equations, v > represents a kinematic viscosity and 
/„ are the Fourier components of the forcing. 
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The three parameters of the model a,b and c are real. In order for the Sabra shell 
model to be a system of the hydrodynamic type we require that in the inviscid (v = 0) and 
unforced (/„ = 0, « = 1,2,3, . . .) case the model will have at least one quadratic invariant. 
Requiring conservation of the energy 

leads to the following relation between the parameters of the model, which we will refer to 
as an energy conservation condition 

a + b + c = Q. (3) 

Moreover, in the inviscid and unforced case the model possesses another quadratic invari- 
ant 

w.|(^)"kP. 

The Sabra shell model ([0 has the following 6 parameters: v,X,ko,a,b, and c. However, 
the "characteristic length-scale" ' does not appear on its own, but only in the following 
combinations: koa,kQb, and hoc. Therefore, without loss of generality we may assume that 
ko = 1. Next, by rescahng the time 

ts-at, 

and using the energy conservation assumption ^ we may set 

fl = 1, /? = — e, c = e — 1. 

Therefore, the Sabra shell model is in fact a three-parameter family of equations with 
parameters v > 0, e, and X > 1 . In most of the numerical investigations of the shell models 
the parameter X was set to A- = 2 (see (31, II14I ). The physically relevant range of parameters 
is \a/c\ > 1, or equivalently, < £ < 2 (see 1141 for details). For < e < 1 the quantity W is 
not sign-definite and therefore it is common to associate it with the helicity - in an analogy 
to the 3D turbulence. The 2D parameters regime corresponds to 1 < £ < 2, for which the 
quantity W becomes positive. In that case the second conserved quadratic quantity W is 
identified with the enstrophy - in analogy to the 3D turbulence. 

Classical theories of turbulence assert that the turbulent flows governed by the Navier- 
Stokes equations have finite number of degrees of freedom (see, e.g., jS), 1121 ). Arguing in 
the same vein one can state that the sabra shell model with non-zero viscosity has finitely 
many degrees of freedom. One of the ways to interpret such a physical statement mathe- 
matically is to assume that the number of degrees of freedom of the model corresponds to 
the Hausdorff or fractal dimension of its global attractor. In our previous study of the Sabra 
shell model of turbulence (|5j|) we proved the existence of a global attractor for the model 
and provided explicit upper bounds of its Hausdorff and fractal dimensions. Therefore, 
we proved that indeed the long-time dynamics of the sabra shell model with non-zero vis- 
cosity has effectively finitely many degrees of freedom. The question remains how many? 
The main motivation behind this work is to provide a lower bound for the Hausdorff and 
fractal dimensions of the global attractor. Namely, to show that for the particular choice 
of the forcing term, and for all £ e (0,2), £ 7^ 1, the Hausdorff and fractal dimensions of 
the global attractor are large, proportional to the upper bound obtained previously in Q. 
However, we also give an example of the forcing such that for £ G (0, 1) and any non-zero 
viscosity v, the long-time dynamics of the Sabra shell model of turbulence is trivial. 

In our work, we show that the Sabra shell model of turbulence possesses a global attrac- 
tor of large dimension for all values of the parameter £ G (0,2), £ ^ 1. In other words we 
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show that for every £ 7^ 1, the Hausdorff and fractal dimensions of the attractor are propor- 
tional to log^^v^' for small enough viscosity v. Therefore, we extend and give a rigorous 
analytical justification for the numerical results observed in 1191 and 1201 for e = 1/2, 
and e = 3/2, which corresponds respectively to the purely "three" and "two-dimensional" 
values of parameters. 

Moreover, in Section |3 we obtain an estimate of the dimension of the global attractor 
in terms of the non-dimensional generalized Grashoff number G defined as 

Ifl 

(4) 



where |f | is an appropriate norm of the forcing term, which will be defined later. More 
specifically, we show that for every e G (0,2), e 7^ 1, and for a small enough viscosity v, 
there exist positive constants ci,C2,c^, depending on A-,e, and independent of the viscosity 
V and the forcing term f, such that 

ci\ogxG + C2 < dim//(.;^) < dimpi^) < ^logxG + c^. (5) 

The right hand side of the inequality was proved in 13, and is true for every forcing term 
f. In this work we show that this estimate is tight in the sense that for particular choices of 
the forcing term f the lower bound in Q is achieved. 

Furthermore, in Section ITTI we study the linear stability of the stationary solution of 
the Sabra shell model, concentrated on a single mode A^. We show that it becomes unstable 
for every and for small enough viscosity for all e G (0,2), e 7^ 1, thus correcting a result 
of (n]. By considering a stationary solution concentrated on an infinite number of shells, 
we are able to demonstrate exactly how the transition to chaos occurs both in the "two" and 
"three-dimensional" parameters regime, through successive bifurcations, as the viscosity v 
tends to zero. 

In the "three-dimensional" regime e G (0, 1), when the parameter e becomes close to 
or 1 the scenario of the transition to chaos is different than in the rest of the interval. 
Namely, for a fixed viscosity, when e crosses the values 0.05 and 0.97, the dimension of 
the unstable manifold of the certain stationary solution drops by a factor of 3. However, the 
attractor in those regimes is still of the size proportional to log^^v^^ the chaotic behavior 
in the vicinity of the particular stationary solution changes dramatically. 

Finally in Section |5] we show that in the "two-dimensional" parameters regime the 
Sabra shell model has a trivial attractor reduced to a single equilibrium solution for any 
value of viscosity v, when the forcing is applied only to the first shell. This result is 
similar to the one for the 2-dimensional NSE due to Yudovich 1221 and independently by 
Marchioro p5l. 

The transition to chaos in the GOY shell model of turbulence was studied previously in 
I3I and 1101 by investigating numerically the stability properties of the special stationary 
solution corresponding to the single mode forcing, which has a '/^ Kolmogorov's scaling 
in the inertial range. It was found that this solution becomes unstable at e = 0.3843, and at 
some value of e the phase transition occurs, when many stable directions become suddenly 
unstable. In this work we show that the nature of the transition to chaos strongly depends 
on the type of the forcing chosen. 

First, we give a brief introduction to the mathematical formulations of the Sabra shell 
model problem. More details on this subject could be found in 13 and |6l. 
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2. Preliminaries and Functional Setting 
In this work we will consider the real form of the Sabra model 

= (afc„+lM„+2M«+l +M„M„+1M„_1 +C^„-lM„_lM„_2) -V^^M„ +/„, 

for « = 1,2,3, .. ., and m„, /„ are real for all n. This formulation is obtained from the 
original one by assuming that both the forcing /„ and the velocity components u„ in the 
equation Q are purely imaginary. Our goal in this work is to show that the upper bounds 
of the Hausdorff and fractal dimensions of the global attractor of the Sabra shell model 
obtained in fS l are optimal in the sense that they can be achieved for some specific choice 
of the forcing. Therefore, this formulation of the model is not restrictive, as long as we are 
able to show in that case that the size of the global attractor matches the upper bound of 

m 

Following the classical treatment of the NSE and Euler equations, and in order to sim- 
plify the notation we write the system Q in the following functional form 

— +vAu + B(u,u) =f (6a) 
dt 

u(0)=u'", (6b) 

in a Hilbert space H. The linear operator A as well as the bilinear operator B will be 
defined below. In our case, the space H will be the sequences space over the field of 
complex numbers M. For every u,v £ H, the scalar product (•, ) and the corresponding 
norm | • | defined as 

(U,V) = l^nVn, |u| = ( ^|m„P 

n=l \n=l / 

The linear operator A : Z)(A) ^ // is a positive definite, diagonal operator defined 
through its action on the sequences u = (mi , M2, . . . ) by 

Au = (k\u\,k^U2, ...), 

were the eigenvalues satisfy the equation (|2}. Furthermore, we will need to define a 
space 

y:=Z)(Ai/2) = {u=(Mi,M2,"3,...) : f;^2|„;f<-}. 

7=1 

The bilinear operator B(u, v) = (Z?i (u, v),Z?2(u, v), . . . ) will be defined formally in the 
following way. Let u = (mi , M2, . . . ), v = (vi , V2, . . . ) be two sequences, then 

£„(u,v) = -k„[Xvn+2u„+i - ev„+iu„^i -X"'M„_iv„_2 + eA."'v„-iM«-2 ' 



for « = 1 , 2, . . . , and where mq = u-i = vq = v_ i = 0. It is easy to see that our definition of 
B(u, v) is consistent with Q. In |5| we showed that indeed our definition of B(u, v) makes 
sense as an element of H, whenever u G H and \ eV oiueV and \ G H. 
For more details on the material of this section see 1 5J and J|6J . 

3. Lower bounds for the dimension of the global attractor - the 
"two-dimensional" parameter regime 

The Hausdorff and fractal dimensions of the global attractor of the evolution equation 
are bounded from below by the dimension of the unstable manifold of every stationary 
solution (see, e.g., Q, 1181 ). Therefore, in order to derive the lower bound for the Haus- 
dorff and fractal dimensions of the global attractor of the Sabra shell model equation we 
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will construct a specific stationary solution of the equation (|6j and count the number of the 
linearly unstable directions of that equilibrium. The same technique was first used in 1161 
(see also [l], I13II ) to obtain lower bounds for the dimension of the Navier-Stokes global 
attractor in 2D. In this section we will consider the "two-dimensional" parameters regime 
of the Sabra shell model corresponding to 

1 < e < 2. 



Consider the forcing 
where 

for 



f=(/l,/2,/3,...), 

r fc", n = mod 3, 
~ \ 0, otherwise, 

1 e-1 5 
a=3log.^ + 3. (8) 

In order to avoid the questions of the existence and uniqueness of the solutions to the 
problem (|6j (see ||5l for details) we will chose some large number M > such that /„ = 0, 
for all « > M. More precisely, such a forcing is supported on the finite number of modes, 
and therefore, according to the results of |5l, for every initial conditions u(0) e H the 
unique solution of the Sabra shell model of turbulence exists globally in time, and possess 
an exponentially decaying dissipation range (see (Sj for details), in particular u{t) e V for 
all f > 0. It was also established in Js] that for such a forcing f the Sabra shell model of 
turbulence has a global attractor, which is a compact subspace of the space V. Later in this 
section we will specify how large the number M should be. 

The corresponding stationary solution of the Sabra shell model equation Q or (|6|l is 

U= (mi,M2,"3,---): (9) 

with 

« = mods, ^^^^ 
otherwise. 

Consider v = (vi, V2, V3, . . . ) e // - an arbitrary perturbation of the stationary solution 
u. Plugging u + V into the equation of motion (|6} we find that the perturbation v satisfies 
the equation 

d\ 

— +vAv + B(u,v)+B(v,u)+B(v,v) =0. 
dt 

To study the linear stabiUty of the equilibrium solution u we will consider the properties 
of the linearized equation 

T 

— +Luv = 0, 
dt 

where the linear operator is defined as 

LuV = vAv + B(u,v)+B(v,u). (11) 

We are looking for the solution of the eigenvalue problem 

LuV = -ov, (12) 

for some o e C. Our goal is to count the number of the solution of equation il2\ with 
Re (a) > 0. The equation il2i in the componentwise form can be written as 

Vk^V,, - k„(Xu„+2Vn+l - eUn+lVn-l + - 1 ) A." ' M„- 1 V„_2 + 

+ A.v„+2M«+i -ev„+iM„-i + (e- l)A,"'v„_iM„_2) = -ov„. (13) 
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where m„ is specified in dlOl i. Note, that u,, = for all n 7^ mod 3, therefore the last 
equation could be written in the following detailed form 

• For n = mod 3, 

vkf,v„ = ~av„. (14) 

• For n = 1 mod 3, 

vfc,^v„-fc„((A.M„+2-eM„-i)v„+i + (e-l)A."'M„_iv„_2) = -ov„. (15) 

• For n = 2 mod 3, 

vk^Vn - k„(Xu„+iv„+2 + {{£ - l)A."'M„-2-eM«+i)v„-i) = -ov„. (16) 
Note that from the relation it follows that o ~ ^^^q correspond to the eigenvectors 

v=(0,...,0,v„„,0,...), (17) 
with v„(, 7^ for some «o = mod 3. However, we are only interested in the solutions of 
the equation Mil with Re{a) > 0. Based on the above the only solution of the relation (I13> 
for which Re (a) > should satisfy v,, = 0, for all n = mod 3. The equations il4l are 
not coupled with the rest of the recursive equations ( I15t and ( I16t . Therefore, in looking 
for non-trivial solutions v of the equation il2\ we can look only for the coupled recursive 
Unear equations il5\ and il6\ . and set 

v„ = 0, V« = mod 3, (18) 

as the solution of ( I14> . 

In what follows we will find sufficient conditions for the existence of non-trivial solu- 
tions for J12t with Re{a) > 0. Denote 



(e- l)^„_iM„_i 



for all « = 1 mod 3, and 



bn.2 



kn+\Un+\ 

vkl+<5 ' 

(e- \)X^^u„-2-^u„+\ 



kfj 7 9, ; 

vkl + a 

for all n ~2 mod 3. Then we can rewrite equations ( I15t and ilbi as a recursive relation 
for Vn 

v„ - b„,iv„+i - c„jv„_2 = 0, for all n = 1 mod 3, 
Vi2—b„,2V„+2^c„2V„-i=0, for all « = 2 mod 3. 

Due to our choice (|8} for the value of a, one can realize from ilQi that c„ 2 = 0, for all 
n = 2 mod 3. Therefore, we can further simplify the last equations, which become 

v„-^„,iv„+i -c„,iv„_2 =0, n = l mod 3, 

Vn - b„.2V„+2 =0, « = 2 mod 3 . 

The following result gives a sufficient condition for the last recursion to have at least one 
non-trivial solution. 
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Lemma 1. Let M be a large positive integer. Let us fix N < M, and assume that N = I 
mod 3. Then the recursive equation M91 has a non-trivial solution of the form v„ ~ 0, for 
all n > A^, and v„ ^ 0, for some n < N, if and only if 

ca^.i = 1. (20) 

Proof. The proof of the Lemma ^ is simple once we observe that the recursive relation 
( I19t for solutions of the form v„ = 0, for n > A^, becomes 

V'l -^I,IV2 = 0, 
V2 -^2,2V4 = 0, 
V4-^4,lV5-C4,iV2 = 0, 
V5 -^5,2V7 = 0, 

VN-2 - bN-2.2VN = 0, 

vn-cnavn-2 =0. 

The last two equations have a one-parameter family of nontrivial solutions if and only if 
the condition i2Q\ is satisfied. □ 

Finally, we are ready to prove the main result of this section. 

Theorem 1. The Hausdorff and fractal dimensions of the global attractor of the equation 
(|6} in the parameter regime 1 < e < 2, with the forcing term f specified in 0, satisfy 

2 1 

dimFH>AimHSA >— ^log;^v-' + — — ^log;,(e-l). (21) 

4-log;,V- 8-21og;^^ 

Proof. Fix M to be large enough, and let < M be such that A' = 1 mod 3. Suppose, that 
for such A^ the condition ((^U) is satisfied for certain o = CJ(A^), depending on A^, for which 
Re (a) > 0. Then, for such o, there exists a solution of equation ( I19> , and in particular, 
there exists a solution of the eigenvalue problem (I12> with Re (a) > 0. Moreover, it is not 
hard to see that if A^i -^N2, then the solutions of the eigenvalue problem (I12> corresponding 
to o(A^i) and o(A^2) are different. 

Therefore, for a given M, in order to count the number of unstable directions of the 
stationary solution ilQ\ . we need to count the number of A'-s, such that (i) N < M; (ii) 
N = I mod 3; (iii) A^ satisfies (I20> with the eigenvalue a, for which Re {a) > 0. 

Let us fix A^ > 0, satisfying N — I mod 3. The condition i20\ becomes 

kN+iUN+i (e— l)^iV+iMA'+i _ 

We get the quadratic equation in a 

0^ + (v4 + v4+2)o + v^44+2 - (e- l)4+i4+i =0- 
This equation has a real positive solution, provided 

v^44+2 - (e - l)4+i4+i < 0. (22) 
Substituting (llOt we obtain the equivalent condition to i22i 

I 1 \ -2i 2(ct-2) 2 1 2 

(e-l)v %Vi >v 4+1. 
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Rearranging terms, the following conditions guarantees the existence of a positive real 
eigenvalue for (I12t 

(8-i)'/V2>x(3-«)(^+i). 

Now, we substitute the value of a from (|8} to obtain 



(iV+l) 



Finally, we get the estimate 



3 



4-log;, V 

6 . 3 



, ,-log^v-' + — — ^log,(e-l). (23) 
4-log;,V 8-21og;,V 
Therefore, we showed that if the M that we have chosen at the beginning of the proof, is 
larger than the right hand-side the relation ( I23> , then for such a choice of the forcing term, 
the number of the unstable direction of the stationary solution dlOt is bounded from below 
by 

2 1 

— — log;, V- ' + — — — log;, (e - 1 ) , 

4-logxV 8-21og;,V 
and the statement of the theorem follows. □ 

In (5) we showed that the dimension of the global attractor of the Sabra shell model 
of turbulence is proportional to the log;,v^' for small enough viscosity v. Therefore, our 
result proves that this bound is tight. 

4. Lower bounds for the dimension of the global attractor - the 
"three-dimensional" parameters regime 

The result obtained in the previous section did not give an answer for the case 

0<£<1, (24) 

which is also known as the "three-dimensional" range of parameters. Therefore, we will 
need to apply different strategy. First, we will consider the linear stability of a stationary 
solution, corresponding to the force acting on a single mode number A^, for some N > 0. 
We will show that for every choice of A^, and for every value of the parameter e G (0,2], 
e 7^ 1, such a stationary solution becomes linearly unstable for sufficiently small viscosity 
V. The stability of a single-mode stationary solution was numerically studied previously in 
111 II . where it was stated that such a solution becomes stable around e = 1. Our rigorous 
proof contradicts this numerical observation. 

Next, we will construct a special type of an equilibrium solution, for which we will be 
able to count the number of unstable directions. The draw-back of this method is that we 
are not able to obtain the exact dependance of the bounds on the parameters of the problem 
e and X, as we succeeded in the "two-dimensional" parameters case. 

4. 1 . On the linear stability of a "single-mode" flow. Let us fix > 1 and consider the 
forcing acting on the single mode of the form 

f'^ = (0,...,0,v^„\0,...), (25) 

where all the components of f'^, except the A^-th, are zero. Consider one particular choice 
of an equilibrium solution corresponding to the above forcing 

u^ = (0,...,0,^^^0,...), (26) 
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which is the analog of the Kolmogorov flow for the Navier-Stokes equations. 

Linearizing the equation Q around the equilibrium solution u'^ and writing the equa- 
tion il2\ in the component form we get the following set of equations. For every j £ N, 
satisfying 2 < \ j —N\, or j — N, we have 

a. 



vk^vj = -avj, 



(27) 



accompanied with the four equations, coming from the nonUnear interaction with u'^ 



_iVN-i ~kNk]^^VN+i+ekN-ik^^vN-2 = -<^vn-i 
+ (1 -e)kNk^^vN-i+ekN+ikf^^vN+2 = -ovn+i, 

V^A'+2^A'+2 + (1 - e)kN+lkf^^\'N+l = -OVn+2- 



Therefore, the eigenvalues of the linear operator L^w (see il l\ ) are a = vk), for 2 < 
\j ~N\, or for j = N, corresponding to the eigenvectors v = (0,0, . . . , 1,0, . . . ), with 1 at 
y-th place. Clearly, those eigenvalues are positive, corresponding to Re{(3) < 0, therefore 
they do not contribute to the number of the linearly unstable directions of the equilibria u^. 
Other eigenvalues of the linear operator L„(v are the eigenvalues of the following matrix 



Jn - 



( v4-2 


% A, 









(l-e)^iv' 


V 








v4+i 






v4+2. 



(28) 



. ) of the linear operator 



(l-e)^;v' 

which will correspond to the eigenvectors v = (vi,V2,V3,. 
with the only non-zero components v^, < \ j —N\ < 2. 

Our goal is to find the condition on the parameters A^, e, and v, such that the matrix Jn 
has eigenvalues with the negative real part, which will correspond to o satisfying Re {a) > 
0. Let us rewrite the expression j28> it in the following way 



JN = k 



-2 



-I 



V 



(1-8) 







-1 

(l-e)?i 





eX 



(29) 



where we denoted for simphcity 

P = v4. 

First, by substituting e = 1, we find that for this value of e the eigenvalues of the matrix Jn 
has always positive real part. Therefore, we conclude that in the case of e = 1 the solution 
u'^ is stable for every and any v. 

For other values of the parameters, we substitute X 
polynomial of the matrix Jn 



2, and write the characteristic 



325 
l6" 



P-^ 



4e^ 
/325 

+ - 



19 



1 



P^ 



4497 
257 
239 



P^ 



197 
1 



Pa 

Fp2 



p2 = 0. 



(30) 



16 16 

Next, by fixing an e, we find the largest p when the real part of the roots of the polynomial 
(13 0> changes its sign. The result of this calculation is shown at Figure 1 . 
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e e 
(a) "Thiee-dimensional" parameter regime. (b) "Two-dimensional" pai'ameter regime. 

Figure 1 . The biffurcation diagram P vs. e. The dashed line indicates 
the appearance of one real negative eigenvalue of the matrix J^, which 
happens only in the regime < e < 1 (a). The solid line shows the 
point at which the real part of two conjugate complex eigenvalues of 
become negative. This bifurcation disappears in the "three-dimensional" 
parameter regime (a) at < e < 0.05 and 0.97 < e < 1. Observe that for 
positive viscosity v the solution u'^ is linearly stable for every at e = 1, 
which can be shown rigorously. 



For every e in the "three-dimensional" subrange of parameters 0.05 < e < 0.97, there are 
two bifurcation points. First, there exists a value of P for which one of the real eigenvalues 
of Jn crosses and becomes negative. Further decreasing p we observe another bifurcation 
at which real part of a pair of complex conjugate eigenvalues becomes negative. Therefore, 
for < £ < 1 there exists a function m(e), such that the matrix 7^? has 3 negative eigenvalues 
for 

0<p = v4<'«(e), (31) 

or in other words, for satisfying 

4A?<log^v-'+log;,m(e), (32) 

Note, that m(e) > 0, at e e (0.05,0.97), and m(e) = otherwise. 

For the range of parameters < £ < 0.05 and 0.97 < e < 1 we observe only one bifur- 
cation point at which one of the real eigenvalues becomes negative. Therefore, we call this 
regime - the "non-chaotic" range of parameters and the reason for that will be explained 
further. 

Finally, for the "two-dimensional" range of parameters 1 < e < 2 the scenario is a little 
different, as we observe only one bifurcation point at which real part of a pair of complex 
conjugate eigenvalues becomes negative. Namely, for 1 < e < 2 there exists a function 
d (e), such that the matrix has 2 negative eigenvalues for 



0<p=v4<'^(e), 



(33) 
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or equivalently, for satisfying 

AN<logy^v-^+logy^e{e). (34) 
In this case we also have, e(e) > for all e G (1,2], and d{l) = 0. 

4.2. Calculating the lower bound of the dimension of the attractor. In the previous 
section we show that the Sabra shell model has at most three unstable direction for the 
"single-mode" forcing. Therefore, we need other type of the force to get a large number of 
unstable direction, and finally, the lower-bound for the dimension of the global attractor, 
which would be close to the upper-bound calculated in (Sj . Let us define the forcing 

g^ffSi Igl^v , ^ (35) 



where f^^ is defined in i25l . Then the stationary solution corresponding to that forcing is 

u^=i»'', (36) 

where u^^ is defined in ( I26l l. Using the results of the previous section on the stability of the 
single-mode stationary solution we conclude that for < e < 1, the number of the unstable 
directions of the solution u^ equals to 3N/5, where satisfies the relation i32\ . On the 
other hand, the number of the unstable directions of the solution u^ for 1 < e < 2 equals to 
2N/5, where satisfies relation ( 134^ 

Recall the definition of the generalized Grashoff number (|4}, which in our case satisfies 

|g| 1 



G = 

Therefore, we can rewrite the bounds i32\ and i32\ in terms of the generalized Grashoff 
number to obtain 

4A^ = log;,v-'+log^/(£) <log;,G + log^/(e), (37) 
where /(e) denotes m(e) or (e). Therefore, we proved the following statement. 

Theorem 2. The Hausdorff and fractal dimensions of the global attractor A of the Sabra 
shell model of turbulence with V > and the forcing defined in \35\ satisfies 

AimF{A)>AimH{A)>K\og^G + \og■^f{^), (38) 

for the positive constant K depending on e satisfying 

( ^, for 0.05<e<0.97, 
K{€) = 1 for l<e<2, (39) 
[ ^, for 0<£<0.05, and Q.91 <E<1, 

and some positive real function /(e), which is only for e = 1. 



The lower bounds for the global attractor, given by the last Theorem do not match 
exactly the upper bounds which were obtained previously in [5 1, namely 

dim// {A)< dimf \ogx G - C(e) , (40) 

where the function C(e) stays positive and bounded for every e <E (0,2). Moreover, the 
constant K in front of the \ogx G term, although can be slightly improved, cannot be brought 
much closer to j to match the upper bound of (I40> . 
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5. Existence of a trivial global attractor for any value of v 

It is well known that the attractor for the 2-dimensional space-periodic Navier-Stokes 
equation with a particular form of the forcing can consists of only one function. This 
well-known example is due to Yudovich 1221 and independently by Marchioro |15.J (for 
the proof see also Q). The same is true for the Sabra shell model for 1 < e < 2, therefore, 
we need to stress that the bounds that we obtained for the dimension of the global attractor 
are valid only for the particular type of forcing that we used in our calculations. 

We mentioned in the introduction that for the 2-dimensional parameters regime the 
inviscid Sabra shell model without forcing conserves the following quantity 

\AM^^j:kty\u„\\ 

for 4y = — logj^^(e — 1). For m > we denote by P„, - the projection onto the first m 
coordinates of the sequence u, and Q,„ = I — P,„. 

Theorem 3. Suppose that the forcing f acts only on the N-th shell for some N > I. Let 
u{t) be the solution of the the equation (|6j in the "two-dimensional" regime of parameters 
1 < e < 2, such that for Y = — | log;^^ (e — 1 ) we have 

Re{B{u,u),A^''u) ^0. 

Then 

limsup|e,„H(f)|2 < c^liminf|P,„H(f)P, (41) 



Ic 



4y 

m+l 



t- 



for C = ^ ■^_4y and m > N. 

Proof. Taking the scalar product of the equation (|6j with u and with A^^u we get two 
equations 



and 



~\u\^+v(Au,u) =Re{fNU*^), 
i||A%|2+v(Au,A%) =/?e(4V^M]^). 



Multiplying the energy equation by k'^J and subtracting it from the last equation we get 

i|(|A%|2-4V|2)+v(|A^+V2„|2_4r||„||2)^0^ (42) 

On the other hand, 

|A^+'/2„|2_4r||„||2^£^2(^4,_^4,)|„^|2> 

«=1 

n=\ 

Plugging the last expression into (I42t yields 

^(|A%|2-4V|2) < -2v4(|A%|2-4V|2), 

and therefore, 

limsup(|A^u(f)p-4V(0l') =0. (43) 
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Finally, 



|Q„,u|2= ^ |«,f = £ ^_^|«,|2< 



< ,4, ^_ 4, (|Q,„A^up-4^|Q,„u|^) = 

%7+l % 



and the result follows from (I43t . 

□ 

Corollary 1. r/ze global attractor of the Sabra shell model of turbulence in the "two- 
dimensional" regime of parameters 1 < £ < 2 with the force applied only to the first shell 

/ = (/,0,0,...), (44) 
is reduced to a single stationary solution 



«' = (-(2,0,0,...). 
vk\ 



|vi(f)|2+v^?|yi(f)P + vi(r)M2(f)«3(f)=0. 



Proof. Let u = (m 1 , M2 , ■ ■ ■ ) be a solution of the Sabra shell model with the forcing f defined 
by J44> . Then it immediately follows from Theoreml^that 

limsup|Qju|^ = 0, 

which means that lim sup,^„ | m„ | = 0, for every n>2. 

Define v = (vi,V2,...) asv = u — u', which satisfies the equation 

d\ 

hvAv + B(u,u) =0, 

dt 

where we used the fact that B(u\u') =0. Taking the inner product of the equation with 
the vector Piv = (vi ,0,0, . . . ) we get that |vi(f)p satisfies 

1 d 

2dt' 

Using the fact that U2{t),Ui{t) tend to as f ^ 00 we conclude that |vi(f)p ^ as f ^ 00. 
Therefore, 

limsup|v|^ = limsupju — u' |^ = 0. 
finishing the proof. □ 

6. Conclusion 

In this work we obtained lower bounds for the dimension of the global attractor of the 
Sabra shell model of turbulence for specific choices of the forcing term. Our main result 
states that for these specific choices of the forcing term the Sabra shell model has a large 
attractor for all values of the governing parameter e € (0,2) \ {1}. We also showed the 
scenario of the transition to chaos in the model, which is slightly different for the two- and 
three-dimensional parameters regime. In addition, in the three-dimensional parameters 
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regime, e e (0, 1 ), we found that when the parameter e becomes sufficiently close to or to 
1 where the chaotic behavior in the vicinity of the stationary solution changes dramatically. 

Finally, we show that in the "two-dimensional" parameters regime the Sabra shell model 
has a trivial attractor reduced to a single equilibrium solution for any value of viscosity v, 
when the forcing is applied only to the first shell. This result is true also for the two- 
dimensional NSE due to Yudovich |2,2I| and independently by Marchioro 1151 (see also 
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